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Grant  AFOSR-83-OI 07 

Interplanetary  Dust  and  the  Visible/Infrared  Sky  Background  Radiation 
Background  Information 

There  is  only  one  source  of  information  on  the  integrated  properties  of 
the  "steady  state"  small  particle  ensemble  of  interplanetary  objects:  the 
zodiacal  light.  Although  brightest  in  the  zodiacal  band  (hence,  its  name), 
the  zodiacal  light  extends  over  the  entire  sky  and  is  caused  by  sunlight 
scattered  and  absorbed  by  Interplanetary  grains.  Observations  of  brightness, 
polarization,  and  color  of  the  scattered  radiation  and  of  the  thermal  emission 
contain  information  on  size,  shape,  composition,  and  overall  properties  of  the 
grains.  Optical  observations  made  from  different  locations  in  the  solar 
system  and  infrared  observations  at  different  wavelengths  from  any  location 
give  the  change  in  zodiacal  light  with  heliocentric  distance.  These  data  are 
needed  to  determine  the  spatial  distribution  of  the  dust  and  changes  in  its 
optical  properties  as  function  of  sun-particle  distance  and  height  above  or 
below  the  plane(s)  of  maximum  concentration. 

In  terms  of  the  scattering  or  photometric  characteristics,  the  zodiacal 
light  cannot  be  seen  "by  itself"  except  close  to  the  sun  where  it  is  very 
bright  compared  to  the  background  starlight.  At  visual  wavelengths  in  and 
near  the  galactic  plane  the  background  starlight  dominates,  whether  observed 
from  space  or  from  the  ground.  For  ground  observations  near  the  horizon, 
airglow  line  emissions  -  and  in  some  wavelength  regions  even  airglow  continuum 
emission  -  dominate.  Infrared  observations  from  the  ground  are  generally 
limited  to  "dry"  sites  and  to  selected  infrared  windows  (spectral  regions 
where  infrared  transparency  of  the  atmosphere  is  high). 

Zodiacal  light  was  expected  to  be  the  limiting  factor  for  astronomical 
observations  from  1-20pm  (Report  of  the  Joint  U.S.,  Netherlands,  and  U.K. 
Scientific  Working  Group  on  the  IR  Astronomy  Satellite  -  IRAS,  and  IR 
observations  such  as  those  of  Price,  et  al.,  Astron.  £. ,  85,  765,  1980).  This 
was  confirmed  by  the  recent  IRAS  observations  (see  the  various  results 
reported  in  A£.  £.  Letters,  278,  March  1,  1984  and  Science,  224,  April  6, 
1984). 
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The  IR  output  of  the  interplanetary  zodiacal  dust  cloud  can  only  be 
understood  in  the  full  context  of  the  physics  of  light  scattering.  Sunlight 
warms  the  circumsolar  dust  to  room  temperature  (near  the  Earth),  causing  it  to 
emit  IR  radiation.  The  rate  of  particle  heating  is  determined  by  the 
particle’s  ability  to  absorb  the  visual  radiation  supplied  by  the  sun.  The 
heating  process  is  a  function  of  particle  size,  shape,  visual  index  of 
refraction  and  its  heliocentric  distance  -  properties  which  can  only  be 
inferred  from  laboratory  studies,  theory,  and  multicolor  observations  from 
1  AU  and  from  a  deep-space  probe.  Some  of  these  data  exist,  but, 
historically,  more  time  and  effort  have  been  devoted  to  morphological  studies 
of  the  zodiacal  light  than  to  analysis  in  terms  of  the  scattering  and 
absorption  of  sunlight  by  interplanetary  dust  grains... and  there  are  still 
large  gaps  even  in  the  morphology: 

-  angular  extent  (most  observations  have  been  made  in  or  near 
the  ecliptic) 

-  color  (observations  have  generally  been  made  in  one  or  two 
"nearby"  wavelengths,  1/A  not  very  different) 

-  resolution  in  sky  coverage  is  limited. 

In  addition,  polarization  and  spectroscopic  observations  are  relatively  rare, 
and  data  for  near  sun  (closer  than  30°)  regions  are  especially  lacking. 

Space  Astronomy  Laboratory  staff  have  studied  the  diffuse  astronomical 
background  radiation  for  over  25  years,  with  emphasis  on  increasing 
observational  accuracy  and  sky  and  wavelength  coverage  and  on  deconvolving  the 
zodiacal  light  to  provide  information  on  the  interplanetary  dust.  Thi3  report 
covers  18  months  (April  1983  through  September  1984)  of  a  continuing  long-term 
study  of  interplanetary  dust  and  the  visible/infrared  sky  background 
radiation. 


Overall  Program  Objectives  and  Stratej 


A  multidisciplinary  approach  is  developed  to  provide  information  on  the 
optical  and  physical  properties  of  interplanetary  dust,  its  location  in  space, 
its  origin  and  dynamics,  and  its  contribution  to  the  astronomical  background 
radiation  in  the  infrared.  This  approach  makes  use  of  a  large  library  of 
observational  data  on  the  brightness,  polarization,  and  color  (near  UV  to  near 
IR)  of  the  zodiacal  light,  including:  extensive  ground-based  observations  of 


the  light  of  the  night  sky  from  Mt.  Haleakala,  Hawaii,  the  first  observations 
of  the  change  in  zodiacal  light  with  heliocentric  distance  (Pioneer  10  and  11 
deep-space  probes),  the  first  multicolor  observations  of  large  regions  of  the 
sky  from  space  (Skylab),  and  additional  sky  and  wavelength  coverage  from 
Shuttle  mission  3  (March  1982).  We  have  reduced  and  analyzed  some  of  the 
Pioneer  10  and  Skylab  data  and  combined  these  data  with  ground-based  results 
to  map  zodiacal  light  brightness,  polarization,  and  color  over  selected 
regions  of  the  sky  as  a  function  of  heliocentric  distance.  These  optical 
(scattering)  data  contain  unique  information  on  the  nature  and  distribution  of 
these  .nterplanetary  grains.  Parallel  studies  have  been  made  of  the  dynamics 
of  small  particles  in  the  solar  system;  of  the  scattering  properties  of 
various  particle  sizes,  shapes,  compositions,  and  spatial  distributions  [the 
Space  Astronomy  Laboratory  has  a  long-term,  AR0/AF0SR  -  supported  program  of 
theoretical  studies  and  microwave  analog  scattering  measurements  to  determine 
the  scattering  by  irregularly  shaped  particles];  and  of  method0  for  inverting 
the  zodiacal  light  brightness  integrals  (scattering  and  emission)  to  derive 
information  on  the  spatial  distribution  and  properties  of  the  grains. 

Our  ultimate  goal  is  to  characterize  the  interplanetary  dust  complex  via 
its  visible  and  IR  3ky  background  radiation  -  especially,  to  evaluate  and 
understand  existing  IR  observations  of  zodiacal  emission  and  to  predict 
brightness,  color,  angular  distribution,  spectral  signatures,  and  structure  in 
observations  not  yet  available. 


Summary  of  Results  (After  18  months) 


-  The  Mt.  Haleakala,  Hawaii  night  sky  observing  facility  and 
photopolarimeter  were  reactivated  and  used  in  an  observational  search  for 
the  Schuerman  Dust  Arcs.  Observing  conditions,  celestial  aspect,  and 
limited  resources  resulted  in  inadequate  data;  the  search  should  be 
repeated. 

-  A  comparison  of  laboratory  microwave  scattering  data  and  photopolarimetry 
of  the  tail  of  Comet  Ikeya-Seki  (1965  VIII)  suggests  that  only  a 
relatively  small  set  of  particle  material/shape  parameters  could  account 
for  the  observed  cometary  polarization. 
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Studies  of  interplanetary  dust  dynamics  included  the  effects  of 
perturbations  by  the  inner  planets  and  the  gravitational  zone  of  influence 
of  a  planet  acting  on  small  celestial  bodies. 

A  new  analytical  inversion  formalism  (Henyey  -  Greenstein)  was  developed 
for  use  with  the  zodiacal  light  brightness  integral  to  derive  the  volume 
scattering  phase  function  of  interplanetary  dust. 

Pioneer  10  and  ground-based  observations  of  zodiacal  light  brightness  were 
found  to  have  small  scale  structures,  with  the  former  structures 
reproducing  at  the  same  elongations  above  and  below  the  ecliptic  and  at 
different  heliocentric  distances. 

A  relatively  high  resolution  observational  model  of  the  zodiacal  light  at 
5080A  was  prepared  for  elongations  30°  to  56°  and  ecliptic  latitutdes 
-30°  to  30°. 

Pioneer  10  observations  of  zodiacal  light  seen  from  beyond  the  Earth 
(between  1  and  3  AU)  were  used  to  depict  changes  in  the  topology  (shape) 
of  the  zodiacal  light  with  sun-spacecraft  distance  and  with  angular 
distance  from  the  sun.  Comparisons  were  made  with  the  corresponding  views 
from  the  Earth  and  from  the  Helios  1  and  2  spacecraft. 

These  results  are  described  in  detail  in  grant-related  publications 
(see  list  which  follows  and  the  attached  copies  of  the  first  page  of  each 
of  these  publications). 
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AN  OBSERVATIONAL  SEARCH  FOR  THE  SCHUERMAN  DUST  ARCS 


F.  Giovane,  J.  L.  Weinberg,  H.  M.  Mann,  and  J.  P.  Oliver 
Space  Astronomy  Laboratory,  University  of  Florida 
l8l0  NW  6th  Street 
Gainesville,  Florida  32609 
U.S.A. 


ABSTRACT.  Large-scale  interplanetary  dust  arcs  were  predicted  by 
Schuerman  (i960)  to  be  associated  with  the  planets.  An  arc,  if  one 
exists,  would  produce  asymmetries  on  opposite  sides  of  the  Gegenscheln 
during  the  Earth's  two  passages  through  the  arc  each  year.  A  brief  des¬ 
cription  is  given  of  Gegenschein/arc  observing  constraints,  of  future 
viewing  opportunities,  and  of  the  results  of  a  test  program  conducted 
from  Mt.  Haleakala,  Hawaii,  in  Jan/Feb  1984  when  the  Earth  encountered 
the  arc  associated  with  Saturn. 

In  the  course  of  performing  dust  dynamics  studies  associated  with  the 
International  Solar  Polar  Mission  (ISPM)  Zodiacal  Light  Experiment 
(Schwehm,  et  al.,  1981)  and  with  analysis  of  Pioneer  10/11  data, 
Schuerman  (19&0)  found  that  a  previously  unsuspected  solar-system 
phenomenon  might  exist:  interplanetary  dust  arcs  which  could  span  the 
entire  solar  system.  An  arc,  hereafter  called  8chueman  Dust  Arc  or  EDA, 
was  predicted  to  be  associated  with  each  planet  -  those  of  Jupiter  and 
Saturn  being  predominant.  His  prediction  was  based  on  the  results  of  re¬ 
calculating  the  restricted  three-body  problem  to  include  radiation  pres¬ 
sure.  The  combined  effects  of  sunlight  and  gravity  extend  the  classical 
equilibrium  (Lagranglan)  points  L4  and  L5  Into  a  circular  arc  extending 
from  L4  through  the  sun  to  -  an  arc  made  up  of  "trapped"  .01  pm  to 
1.0  Mm  size  particles.  For  particles  in  this  size  range,  the  radiation 
pressure  force  is  not  negligible  compared  to  the  gravitational  force. 
Such  grains  are  referred  to  as  6  particles.  Perturbing  and/or  destructive 
forces  such  as  the  Lorentz  force  resulting  from  the  interplanetary  mag¬ 
netic  field,  solar  wind  interaction,  and  gravitational  effects  of  other 
planets  were  not  included  in  Schuerman1 s  treatment.  Thus,  the  stability 
of  the  arcs  has  not  been  definitively  established.  Schuerman  estimated 
that  for  a  typical  value  of  6  ■  0.57,  particles  in  the  Juplter-sun  and 
flaturn-sun  arc  systems  have  characteristic  lifetimes  T  of  14,700  and 
$0,000  years,  respectively.  Since  8  particles  now  exist  in  the  inter¬ 
planetary  median  and  the  gradual  aplrallng  into  the  sun  due  to  the 
Poyntlng-Robertaon  effect  occurs  on  a  time  scale  on  the  order  of  T,  then 
the  replenishment  mechanism,  regardless  of  details  as  to  its  source, 

muet  be  of  sufficiently  short  time  scale  to  replenish  the  interplanetary 

is 
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LABORATORY  RESULTS  ON  POLARIZATION  PROPERTIES  OF  ELONGATED  PARTICLES 
AND  COMPARISONS  TO  DUST  IN  THE  TAIL  OF  COMET  IKEYA-SEKI  (1965  VIII) 


BO  X.  S.  GUSTAFSON 

Space  Astronomy  Laboratory,  University  of  Florida 
1810  NW  6th  Street 
Gainesville,  Florida  32609 
U.S.A. 


ABSTRACT.  Data  obtained  through  microwave  analog  measurements  are  used 
to  compare  optical  polarization  characteristics  of  silicate  particles 
ranging  in  shape  from  spheroids  and  cylinders  to  chains  of  spheres*  It 
is  confirmed  that  the  dependence  of  the' degree  of  polarization  on  scat¬ 
tering  angle  Is  a  powerful  indicator  of  particle  shape.  The  reversal 
and  steep  gradient  in  the  degree  of  polarization  found  by  Weinberg  and 
Beeson  (1976a, b)  in  the  dust  tall  of  comet  Ikeya-Seki  (1965  VIII)  is 
suggestive  of  laboratory  measurements  corresponding  to  1.6  pm  long, 

0.4  pm  in  diameter  silicate  cylinders. 


1 .  INTRODUCTION 

Krishna  Swamy  (1978)  showed  that  the  run  of  the  degree  of  polarization 
with  scattering  angle  observed  in  the  tail  of  comet  Ikeya-Seki  (1965 
VIII)  could  result  from  size  segregation  among  spherical  particles. 
Segregation  along  the  tail  is  expected  due  to  faster  diffusion  away  from 
the  tall  for  low  mass  particles.  But  spherical  particles  are  probably 
rare  in  a  disintegration  event.  Cubes  or  other  rough  shapes  may  be  a 
better  approximation.  Sekanina  and  Farrell  (1980)  found  Independent 
evidence  for  elongated  grains  in  the  call  of  comet  West  (1976  VI),  and 
suggested  that  chain  like  aggregates  may  have  been  common.  Schuerman 
ec  al.  (1981)  and  Zerull  (1985)  concluded  that  there  are  marked  devia¬ 
tions  from  Mle  polarization  patterns  for  the  caa#  of  non-spheres.  Tha 
degrees  of  polarization  produced  by  chains  of  five  spheres,  similar 
sized  circular  cylinders  and  spheroids  are  compared  to  illustrate  tha 
effect  of  detailed  shape  including  sharp  edges.  In  this  paper,  we  con¬ 
centrate  on  the  degree  of  polarization.  This  quantity,  a  ratio  of  in¬ 
tensities,  is  independent  of  the  (unknown)  number  of  particles  along  tha 
line  of  sight.  In  addition,  the  angular  dependence  of  the  degree  of 
polarization  is  found  to  ba  a  powerful  indicator  of  particle  parameters 
(see  also  Schuerman  at  al.,  1981  and  Zerull,  1985). 
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A  NEW  APPROACH  TO  EVALUATE  PLANETARY  PERTURBATIONS  ON  A  CLOUD  OP  DUST 
IN  LOW  ECCENTRICITY  HELIOCENTRIC  ORBITS 


BO  X.  S.  GUSTAFSON 

Space  Astronomy  Laboratory,  University  of  Florida 
1810  NW  6th  Street 
Gainesville,  Florida  32609 
U.S.A, 

ABSTRACT.  Dynamical  perturbations  on  enaanblea  of  particles  in  helio¬ 
centric  orbits  of  low  eccentricity  are  integrated  over  time.  The  dust 
is  perturbed  by  radiation  pressure,  Poyntlng-Robertsoa  drag,  their  cor¬ 
puscular  counterparts,  and  by  gravitation  due  to  any  number  of  planets. 

A  dust  cloud  is  represented  by  a  set  of  centroids  and  orbital  disper¬ 
sions  (about  the  centroids).  Gravitational  perturbations  on  the  cen¬ 
troid  are  derived  from  a  single  matrix,  valid  for  any  planet,  in  the 
appropriate  frame  of  reference.  After  transformation  of  the  time  deriva¬ 
tives  to  a  common  coordinate  system,  the  perturbation  rates  are  summed  up 
and  integrated.  The  time  dependence  of  the  planets'  orbital  elements  art 
evaluated  inside  the  time  integral. 


1 ,  INTRODUCTION  ' 

It  is  time-consuming  and  costly  to  estimate  overall  perturbations  on  s 
dust  cloud  by  computing  trajectories  of  individual  dust  particles.  The 
efficiency  is  increased  through  the  use  of  centroids  to  represent  sets  of 
orbits.  Perturbations  due  to  each  planet  are  computed  separately.  Geo¬ 
metries  are  chosen  such  that  the  planet's  perturbations  on  a  centroid  may 
be  expressed  by  a  single  quantity  depending  only  on  the  angle  that  the 
plane  containing  the  centroid  makes  to  the  planet's  orbital  plane  end  the 
ratio  of  heliocentric  distances.  A  comparatively  simple  interpolation- 
formula  allows  the  matrix  representing  the  perturbations  to  be  stored  on¬ 
line.  Properly  scaled,  a  single  matrix  may  represent  perturbations  by 
any  planet  in  the  appropriate  geometry. 

2.  COMPUTATIONAL  PROCEDURE 
2.1  The  approach 

Perturbations  on  enaemblea  of  particles  are  evaluated  as  shown  in  the 
flow  diagram.  Figure  1.  The  ensemble  is  divided  into  a  set  of  ficti¬ 
tious  orbits,  each  s  centroid  for  a  set  of  duat  trajectories.  Each 
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PLANETARY  PERTURBATIONS:  EFFECTS  ON  THE  SHAPE  OF  A  CLOUD  OF  DUST 
IN  CIRCULAR  HELIOCENTRIC  ORBITS 


BO  X.  S.  GUSTAFSON 

Space  Astronomy  Laboratory,  University  of  Florida 
1810  NU  6th  Street 
Cainesville,  Florida  32609 
U.S.A. 


ABSTRACT.  Synthetic  Infrared  pictures  are  used  to  illustrate  changes  In 
the  shape  of  portions  of  an  interplanetary  dust-cloud.  The  dust  parti¬ 
cles,  in  circular  heliocentric  orbits,  are  perturbed  by  radiation  and 
corpuscular  forces  combined  with  gravitational  disturbances  by  the  major 
planets.  Dust  in  the  inner  solar  system  and  close  to  the  ecliptic  or 
the  orbital  plane  of  Venus  is  brought  to  a  more  narrow  range  of  ecliptic 
lacitudes.  A  dust-band  evolves  near  and  Inside  the  orbit  of  Venus.  The 
cloud's  shape  is  leas  affected  at  high  ecliptic  latitudes. 

1 .  INTRODUCTION 

Derivations  of  the  shape  of  the  zodiacal  cloud  based  on  dynamical  evolu¬ 
tion  are  initial-value  problems.  Because  the  source  remains  unknown  as 
well  as  the  cloud' b  long-term  steady  state  or  transient  nature,  no 
statements  can  be  made  regarding  the  absolute  shape  of  the  interplane¬ 
tary  dust  cloud.  Instead,  we  address  the  problem  of  ohangoo  in  the 
shape  of  a  portion  of  the  cloud  as  the  dust  spirals  through  the  inner 
solar  system.  Conventional  wisdom  states  that  dust  migrates  through  the 
solar  system  while  losing  angular  momentum  to  drag  forces.  It  is  not 
yet  established  whether  transfer  of  momentum  from  the  planets  may  affect 
the  rate  significantly.  The  resulting  uncertainty  in  magnitude  of  the 
perturbations  has  little  effect  on  our  qualitative  results.  Because  or¬ 
bital  eccentricities  are  neglected,  momentum  transfer  and,  by  implica¬ 
tion,  dynamical  life-times  of  the  dust  are  not  addressed  in  this  paper. 

2.  OBSERVATIONS 

Observations  show  that  the  symmetry  plane  of  the  inner  zodiacal  cloud  is 
close  to  the  orbital  plane  of  Venus  (Misconi,  1980,  and  references 
therein).  Leinert  at  al.  (1980)  concluded  that  this  symmetry  extends  to 
1  A.U.,  whereas  most  observations  in  this  region  favor  a  plane  close  to 
the  invariable  plane  of  the  solar  system  (Misconi,  I960).  The  bright¬ 
ness  plane  of  symmetry  is  close  to  the  plane  of  maximum  particle  number 
densities  except  where  shifted  by  the  duet  particles'  light  scattering 
characteristics,  as  in  the  Gegenscheln  (Misconi,  1981). 
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CAN  C0META8Y  DUST  PERTURBED  BY  THE  INNER  PLANETS  BE  AN  EXPLANATION 
FOR  THE  OBSERVED  DISTRIBUTION  OF  INTERPLANETARY  DUST? 

Bo  R.  8.  Gustafson  and  Nebll  Y.  Nisconi 

Space  Aatronomy  Laboratory,  University  of  Florida, 
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ABSTRACT 

Numerical  Integrations  of  perturbation  equations  using  the  Encka 
method  show  that  the  long-tens  dynamical  evolution  of  interplane¬ 
tary  duet  is  strongly  affected  by  the  inner  planets.  Porturbations 
in  inclination  and  ascending  node  nay  explain  the  observed  inclin¬ 
ation  of  the  plane  of  maximum  dust  density.  Non  monotonlc  changes 
ere  found  in  the  semi-major  axes  and  eccentricity.  As  a  result, 
temporal  increases  in  semi-major  axis  aiay  slow  the  sink  of  rod la¬ 
cs  1  dust  into  the  sun.  with  respect  to  current  eatiaetee. 


INTRODUCTION 

Recent  studies  have  shown  that  the  plane  of  maximum  dust  den¬ 
sity  (symmetry  plane)  of  the  zodiacal  cloud  varies  in  inclination 
with  respect  to  the  ecliptic  piano  as  s  function  of  heliocontric 
distance  (Nisconi  and  Weinberg,  1976 i  Hleconl,  1979)  Lelnert,  et 
el.,  1980).  This  suggests  that  we  are  dealing  with  a  multiplicity 
of  symmetry  "planes"  which  tend  towards  the  orbital  planes  of  the 
planets  (Fig.  1).  Consequently,  the  suggestion  was  made  that  long- 
tens  gravitational  perturbations  by  Jupiter  and  the  inner  planets 
are  influencing  the  spatial  density  distribution  of  interplanetary 
dust  (Nisconi,  1977 i  Nisconi  and  Weinberg,  1971). 

Long-torsi  gravitational  perturbations  by  the  inner  planets  on 
the  orbital  elements  of  the  dust  have  not  previously  been  investi¬ 
gated.  Studies  on  the  gravitational  effect  of  Jupiter  (Brouwer 
end  Clemance,  1961 i  end  Smith,  1964)  led  to  the  belief  that  Jupi¬ 
ter  is  the  sole  perturber  of  the  dust  and  that  the  symmetry  plane 
of  interplanetary  dust  is  expected  to  be  near  the  invariable  plane 
of  the  solar  system  which,  in  turn,  is  near  the  orbital  plane  of 
Jupiter.  Since  that  belief  ia  not  supported  by  observations,  the 
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Whereat  the  timer  pianeit'  perturbations  on  meteoroids'  and  larper  interplanetary  bodies'  orbits 
have  been  studied  extensively,  they  are  usually  itefkcied  in  studies  of  the  dynamics  of  smaller 
panicles  producinf  the  aodiacal  lifht  through  scattering  of  sunlight.  Forces  acting  on  these  dusi 
panicles  arc  fairiy  well  known  and  include  radiation  lorccs  and  interaction  with  the  solar  wind.  This 
aniclc  is  the  first  in  a  scries  aimed  at  improving  our  knowledge  of  the  dynamical  evolution  of  dust  in 
interplanetary  space  by  studying  the  combined  effects  of  these  perturbations  including  gravitational 
penurbations  by  the  planets  Venus.  Earth.  Mars,  and  Jupiter.  The  necessity  of  including  effects  of 
the  inner  planets  in  dust  dynamics  investigations  is  established.  Sample  trajectories  are  presented 
lo  iHustrate  commonly  occurring  phe nonsense,  such  as  nonmonotonic  changes  in  semimajor  axis, 
eccentricity,  inclination,  and  in  the  line  of  nodes.  These  perturbations  are  shown  lo  be  due  to  the 
inner  planets  as  opposed  lo  Jupiter  or  nongravitatiotul  forces,  t  tass  AiUrw.  *»„.  tm 


Studies  of  the  zodiacal  light  in  the  last  “planes**  which  tend  toward  the  orbital 
decade  have  supplied  the  first  detailed  in-  planes  of  the  planets  (Fig.  I).  Conse- 
formation  on  the  large-scale  properties  of  quently,  the  suggestion  was  made  that  long- 
the  interplanetary  dust  complex,  and  are  term  gravitational  perturbations  by  Jupiter 
suggestive  of  planetary  influences.  We  and  the  inner  planets  are  influencing  the 
adopt  common  practice  and  refer  to  parti-  spatial  density  distribution  of  intcrplanc- 
cles  with  radii  less  than  100  gun  as  dust,  tary  dust  (Misconi.  1977;  Misconi  and 
This  is  convenient  because  the  large-scale  Weinberg.  1978).  Morfill  and  Grtln  (1979) 
distribution  of  this  size  particles  can  be  have  suggested  that  the  Lorentz  force  re- 
studied  from  observations  of  the  zodiacal  suiting  from  the  interaction  of  the  dust  with 
light,  which  is  produced  through  scattering  the  interplanetary  magnetic  field  is  influenc- 
of  sunlight  by  primarily  10-  to  lOO-jim-size  ing  the  dust  density  distribution.  This  influ- 
dust  (Giese  and  Griln,  1976;  ROser  and  ence  would  give  a  symmetry  toward  the  so- 
Staude.  I978j.  The  spatial  density  distribu-  lar  equatorial  plane.  One  difficult)  with 
tion  and  other  properties  of  the  dust  cloud  Morfill  and  Griln's  suggestion  is  that  the 
were  determined  from  measurements  by  Lorentz  force  becomes  negligible  for  large 
optical,  impact,  and  other  dust  detection  in-  particles  (>  10  pm)  from  which  most  of  the 
struments  on  board  the  Pioneer  10  and  II.  zodiacal  light  is  now  thought  to  arise.  Addi- 
and  the  Helios  A  and  B  space  probes  tional  difficulties  arise  from  uncertainties  in 
(Weinberg  and  Sparrow,  1978).  The  plane  the  value  for  the  charge  of  these  panicles 
of  maximum  dust  density  (symmetry  plane)  and  in  the  nature  and  behavior  of  the  inter- 
of  the  zodiacal  cloud  was  found  to  vary  in  planetary  magnetic  field  through  lifetimes 
inclination  with  respect  to  the  ecliptic  plane  of  the  panicles  of  (he  order  of  HP-10* 
as  a  function  of  heliocentric  distance  (Mis-  years.  Finally,  if  the  Lorentz  force  is  re¬ 
coni  and  Weinberg,  1978;  Misconi,  1979;  sponsible  for  the  observed  symmetry  plane. 
Leinen  et  al,  1980).  This  suggests  that  we  then  one  would  expect  the  plane  to  have  a 
are  dealing  with  a  multiplicity  of  symmetry  smooth  decrease  in  inclination  away  from 
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Summary.  To  determine  the  scattering  characteristics  of  inter¬ 
planetary  particles,  we  substitute  a  linear  combination  of  three 
Henyey-Greenstein  functions  for  the  mean  volume  scattering 
phase  function  in  the  zodiacal  light  brightness  integral.  Residuals 
are  obtained  by  comparing  the  observed  zodiacal  light  in  the 
ecliptic  with  results  of  the  integral.  Minimization  of  the  residuals 
leads  to  a  scattering  function  which  has  a  strong  peak  in  the 
forward  direction,  an  isotropic  part  at  intermediate  scattering 
angles  and  a  mild  enhancement  in  the  backward  direction. 

The  same  method  is  employed  to  analyze  the  polarized 
components  of  the  zodiacal  light.  The  scattering  phase  function, 
*±(0),  for  the  component  polarized  perpendicularly  to  the 
scattering  plane  remains  substantially  higher  than  the  parallel 
component,  over  most  scattering  angles.  But  at  scattering 
angles  greater  than &=*  165°,  Q±(&)  becomes  less  than  4u(d)  with 
a  masimum  difference  between  the  two  occurring  at  05x175*. 

Key  words;  zodiacal  light  -  scattering  phase  function  -  interplane¬ 
tary  dust 
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Understanding  how  interplanetary  dust  particles  scatter  sunlight 
into  different  directions  is  one  of  the  prime  objectives  in  zodiacal 
light  studies  (see  reviews  by  Lcinert,  1975;  Weinberg  and  Spar¬ 
row,  1978).  Such  scattering  characteristics  give  important  clues  for 
the  chemical  and  physical  identities  of  the  interplanetary  particles 
(Giese  el  aL,  1978;  Greenberg  and  Gustafson,  1981 ;  Weiss-Wrana, 
1983).  and  also  provide  an  essential  stepping-stone  for  determining 
the  morphology  of  the  zodiacal  dust  cloud  (Dumont,  1976a; 
Leincrt  cl  aL,  1976;  Buitrago  et  al.,  1981, 1983).  During  the  last 
decode  tenons  attempts  have  been  made  to  determine  the 
taiicnng  phase  function  for  interplanetary  particles  from  the 
observed  brightness  of  zodiacal  light  (Dumont,  1974,  1976a; 
Dumont  and  Sincbcz,  1975b;  Leinerl  et  al,  1976).  Two  principal 
assumptions  employed  in  these  attempts  are  as  follows:  (i)  Scatter¬ 
ing  properties  of  interplanetary  particles  are  independent  of 
location  within  the  zodiaul  dust  cloud,  (ii)  Dust  distribution  in  the 
ecliptic  plane  follows  a  simple  power-law  relation  with  heliocentric 
distance.  There  are  indications  that  these  assumptions  are  unwar¬ 
ranted  (Schuennan,  1980;  Hanner  et  al,  1976);  nevertheless, 
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they  are  necessary  simplifications  taken  in  order  u>  make  the 
problem  manageable.  Thus,  we  should  view  the  resulting  cmptnal 
scattering  phase  function  only  as  an  average  over  rather  impre¬ 
cisely  defined  parameters. 

Most  theoretical  methods  for  denying  the  cm pineal  icaucnag 
phase  function  require  taking  denvatives  of  the  observed  zodiacal 
light  brightness.  As  an  alternative,  we  lake  an  integral  approach  to 
the  problem.  We  utilize  the  versatile  Hcnycy-Gracaucra  funnion 
which  conveys  essential  characteristics  of  a  com  pics 
pattern  by  varying  a  single  parameter.  Substituting  a  linear 
combination  of  Henyey-Greenstein  functions  lor  the  scaucnng 
phase  function  in  the  zodiacal  light  brightness  integral  wc  form 
residuals  between  results  of  the  brightness  integral  and  the 
observed  brightness.  The  residuals  are  then  by  the 

method  of  non-linear  least  squares  to  determine  the  best  combin¬ 
ation  of  parameters  for  the  scattering  phase  function  of  the 
interplanetary  particles. 

In  the  second  section  wc  present  basic  (onnulauom  lor 
deriving  an  empirical  scattering  function,  rmphasinng  the  dif¬ 
ference  betvyfcn  differential  methods  of  inversion  and  our  method 
of  non-linear  least  squares.  In  the  third  section  the  observed  total 
brightness  of  the  zodiacal  light  in  the  cdipuc  is  analyzed  to  show 
advantages  as  well  as  limitations  of  the  integral  method  la  the 
fourth  section  the  polarized  brightness  of  the  zodiacal  light  u 
analyzed  to  determine  polarization  characteristics  of  the  inter¬ 
planetary  particles.  The  paper  concludes  with  ilisnitiiont  oa 
future  improvements. 


L  Inversion  of  the  brightness  integral 

2.1.  Zodiacal  light  brightness  integral 

The  surface  brightness  of  zodiacal  light  Z(s)  seen  in  the  ecliptic  al 
solar  elongation  angle  t  is  the  sum  of  scattered  sunlight  received 
from  all  tfrf  particles  along  the  given  dirocuop,  bests,  Z(<) »  *>*• 
by 

«(»)-  J 1  Flr)air)/M*)«6;i)dtM.  (I) 

B  9 

where  variables  r,  di.  and  8  represent,  as  schematically  ilhiursied 
in  Fig.  1,  the  heliocentric  distance,  line-nkmeni.  and  acniisnag 
angle,  respectively,  and  s  is  particle  radius.  The  talar  ftus  at 
heliocentric  distance  r  is  denoted  by  F(r),  which  can  he  replaced  by 
an  inverse-square  relation  F^rjr)1  with  f,  bong  the  tolar  Out  at 
heliocentric  distance  r0.  Note  that  this  formulation  assumes  the 
particle  scaucnng  properties  to  be  independent  of  location  With 
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ABSTRACT.  A  linear  combination  of  3  Henyey-Greenatein  phase  functions 
la  substituted  for  the  mean  volume  scattering  phase  function  In  the 
zodiacal  light  brightness  Integral.  Results  of  the  integral  are  then 
compared  with  the  observed  brightness  to  form  residuals.  Minimisation 
of  the  residuals  provides  us  with  the  best  combination  of  Henyey- 
Creenstein  functions  for  the  scattering  phase  function  of  zodiacal 
dust  particles. 


1.  INTRODUCTION 

Most  theoretical  methods  for  deriving  an  empirical  scattering  function 
for  Interplanetary  dust  particles  require  taking  derivatives  of  the 
observed  zodiacal  light  brightness.  Observational  errors  are  easily 
amplified  in  the  process  of  differentiation  and  likely  to  distort  the 
scattering  function  derived  therefrom.  As  an  alternative,  we  taka  an 
Integral  approach:  A  utilization  of  the  versatile  Henyey-Creenstein 
(HC)  function  as  a  parameterized  trial  scattering  function  will  enable 
us  to  employ  the  method  of  non-linear  least  squares  in  determining  the 
mean  volume  scattering  phase  function  for  zodiacal  dust  particles. 


2.  BRIGHTNESS  INTEGRAL  WITH  THE  HENYEY-GREENSTEIN  FUNCTION 


A  linear  combination  of  three  HG  functions 


a  -  *»> 
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having  different  asymmetry  factors,  g.  ,  adequately  describes  such 
general  features  as  the  forward  peak,  isotropic  middle,  and  backward 
enhancement  in  a  particle  scattering  pattern.  Here,  w.  represents 
relative  weights  of  the  three  components,  and  8  denotes  scattering 
angle.  With  the  usual  power-law  relation,  n(r)  “  n  (r  /r)  ,  for  tha 
heliocentric  density  distribution  of  particles*  the  zoalacal  light 
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ABSTRACT.  Cegenschein  observations  fro*  Pioneer  10  were  found  to  have 
brightness  structures  with  an  amplitude  of  about  lOt  and  a  period  of 
several  to  ten  degrees  in  elongation.  A  search  is  made  for  such 
structures  in  high  angular  resolution  ground-based  observations  from 
Mt.  Haleakala,  Hawaii.  A  new  empirical  method  la  used  to  correct  for 
atmosphere-originated  radiation.  Background  starlight  is  subtracted 
using  Pioneer  10  observations  from  beyond  the  asteroid  belt.  Prelimina¬ 
ry  analysis  of  the  ground  data  also  indicates  the  presence  of  small 
amplitude  structures  in  the  brightness  distribution. 


This  is  a  first  report  on  our  search  for  small  scale  structures  in  the 
zodiacal  light  brightness.  By  structures  we  mean  small  amplitude  spati¬ 
al  fluctuations  in  the  brightness  distribution.  Such  structures  might 
be  expected  from  inhomogeneities  in  the  spatial  distribution  of  the 
zodiacal  dust  and/or  structure  in  the  scattering  phase  function.  There 
have  been  a  number  of  observational  hints'  for  the  existence  of  such 
structures.  Ground-based  observations  of  zodiacal  light  brightness  by 
Weinberg  (1963)  show  an  enhancement  at  elongation  135*,  which  also  has 
a  counterpart  in  the  polarized  brightness  (Weinberg,  19c*1*;  Frey,  ct  el., 
197 h ) .  Other  space  and  ground-based  observations  of  zodiacal  light  have 
also  shown  irregularities  in  the  brightness  distribution  (Sparrow  and 
Ney,  1972;  Dumont  and  Sdnchez,  1975).  Despite  these  observational  indi¬ 
cations,  structures  in  the  zodiacal  light  have  generally  been  ignored 
for  two  reasons.  First,  there  was  a  lack  of  detailed  information  on  the 
distribution  of  background  sturlight  over  the  sky.  Second,  for  ground 
obaervations  it  is  difficult  to  determine  the  airglow  continuum  an<*  to 
correct  for  airglow  and  astronomical  diffuse  background  that  is  scatter¬ 
ed  into  the  telescope's  field  of  view  by  the  Earth's  atmosphere. 

Figure  1  shows  Pioneer  10  near-Earth  observations  »*f  zodiacal  lignt 
brightness  as  a  function  of  differential  ecliptic  longitude  A-l*  for 
various  ecliptic  latitudes.  Connection  of  the  points  in  Figure  1 
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ABSTRACT 

Tlsserand's  definition  of  the  "sphere  of  action*  of  a  planet  Is. based  on  the  equality  of 
tidal  vs.  gravitational  acceleration  ratios  of  the  sun  and  planet.  Opik  and  others  based 
their  relation  on  equating  the  differential  solar  and  planetary  forces  on  the  particle. 

Neither  espreaaion  was  formulated  to  describe  the  tone  of  Influence  surrounding  a  planet  when 
considering  the  snail  but  significant  (l.e.  long-term)  perturbative  effects  of  the  planets  on 
a  particle's  orbital  elements.  For  the  purpose  of  determining  these  effects  on  interplanetary 
gust  we  dervlve  a  zone  of  Influence  baaed  on  equating  the  gravitational  forces  of  the  sun  and 
planet,  and  demonstrate  Its  applicability  by  utilizing  the  partlole'g  closest  approach  to  the 
planet  as  a  measure  or  the  zone  of  influence, 

INTRODUCTION 

In  the  process  of  computing  the  orbital  evolution  of  dust  partioles  In  Interplanetary 
speee.  It  was  necessary  to  determine  whether  or  not  a  particle  was  close  enough  to  a  planet 
that  integrating  the  perturbing  forces  using  Encke's  method  uss  no  longer  valid.  Both  expree- 
sions  for  the  "sphere"  of  Influence  of  a  planet*  proved  to  be  Inadequate  when  the  geometries  of 
subsequent  Interactions  were  considered.  In  particular,  when  each  relation  was  scaled  to  one 
planet.  It  could  not  be  generalized  to  the  other  planets.  Consequently,  the  formula  for  the 
zone  of  Influence  must  depend  strongly  on  the  specific  subject  to  be  studied.  We  have  studied 
this  dependence  using  chonges  in  several  orbital  elements  as  our  criteria, 

THEORY 

The  general  espresslon  for  the  radial  distance  to  the  boundary  of  the  planetary  zone  of 
Influence  on  small  particles  can  be  written  as 

•  •  ■‘••*p)*g(mp,  other  planetary  parameters)  0) 

where  k  is  a  constant  which  depends  on  the  nature  of  the  problem  under  study,  f  and  g  are 
functions,  a_  i«  the  planet's  semlmajor  axis,  and  mp  is  the  planet's  mass.  The  other  para¬ 
meters  may  include  the  eccentricity  of  the  planet,  as  well  as  the  other  orbital  parameters  and 
such  factors  as  oblateness  and  axial  tut.  Functions  f  and  g  are  generally  not  separable  but 
we  have  found  that  the  erfect  of  semlmajor  axis  can  be  separated  from  the  other  parameters. 

The  equations  of  Tlsserand  0  889)  and  Opik  <1951 )  can  be  written  as 

S  •  (l/M#)2/5epmp2/5.  end  (2) 

a  •V2l'/2Ha)  1/3  *p»p’/3.  respectively.  (3) 

If  we  equate  the  gravitational  and  radiation  forces  of  the  sun  and  planet,  a  different 
relationship  la  obtained.  Consider  a  particle  of  mass  s  at  a  distance  s  from  a  planet  of  mass 
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ABSTRACT .  Tisaerand's  definition  of  the  "sphere  of  action"  of  a  planet 
is  based  on  the  equality  of  tidal  vs.  gravitational  acceleration  ratios 
of  the  sun  and  planet.  Opik  and  others  based  their  relation  on  equating 
the  differential  solar  and  planetary  forces  on  a  particle.  Neither  ex¬ 
pression  was  formulated  to  describe  the  zone  of  Influence  surrounding  a 
planet  when  considering  the  small,  but  significant,  long-term  pertur¬ 
bative  effects  of  the  planets  on  a  particle's  orbital  elements.  For  the 
purpose  of  determining  these  effects  on  interplanetary  dust  we  derive  a 
gone  of  Influence  based  on  equating  the  gravitational  forces  of  tba  sun 
and  planet. 

The  general  expression  for  Che  radial  distance  to  the  boundary  of 
Che  planetary  zone  of  influence  on  small  particles  can  be  written  aa 

a  -  k*f (ap) 'g(mp,  other  planetary  parameters)  (1) 

where  k  is  a  constant  which  depends  on  the  nature  of  the  problem  under 
study,  f  and  g  are  functions  which  may  also  have  the  sane  dependence  aa 
k,  ap  is  the  planet's  semimajor  axis,  and  np  ie  the  planet's  mass.  The 
other  parameters  may  Include  the  eccentricity  of  the  planet,  as  well  as 
other  orbital  parameters  and  such  factors  as  oblateness  and  axial  tilt. 
Functions  f  and  g  are  generally  not  separable  but  we  have  found  chat  the 
effect  of  semimajor  axis  can  be  separated  from  the  other  parameters. 

The  equation  of  TiBserand  (1889)  and  Opik  (1951)  can  be  written  as 

a  -  <1/M8)  '5apmp  and  (2) 

s  »  *j(l/2M„)  ^3a  m  ^3,  respectively.  (3) 
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If  we  equate  the  gravitational  and  radiation  forces  of  the  sun 
and  planet,  a  different  relationship  is  obtained.  Consider  a  particle  of 
mass  m  at  a  distance  s  from  a  planet  of  mass  mp  and  semimajor  axis  a_, 
and  at  a  distance  r  from  the  sun.  The  magnitude  of  the  force  on  the 
particle  due  to  the  planet  la 
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ABSTRACT.  Cround-based  observations  of  the  evening  zodiacal  light  taken 
by  Weinberg  and  Mann  from  Mt.  Haleakala,  Hawaii,  during  March  1966  are 
used  to  derive  a  table  of  zodiacal  light  brightnesses  at  high  spatial 
resolution  (as  little  as  0.5°  in  differential  ecliptic  longitude  1-1 
and  1.0*  in  ecliptic  latitude  8)  over  the  region  29. Sv  1-1  < 56*.  * 

-30*<  0  <  30° .  Significant  differences  are  found  in  the  brightness  dis¬ 
tributions  above  a-  d  below  the  ecliptic  plane. 

Brightness,  polarization,  color  and  angular  dependence  of  the  light 
of  the  night  sky  were  systematically  observed  by  Weinberg  and  Mann  from 
Mt.  Haleakala,  Hawaii  between  1965  and  1969  (Weinberg  and  Mann,  1967). 

One  of  the  several  observing  techniques  that  was  used  involved  scanning 
Che  main  cone  of  the  zodiacal  light  over  a  range  of  160*  in  azimuth  cen¬ 
tered  on  the  ecliptic,  beginning  (evening)  or  ending  (morning)  with  the 
onset  of  astronomical  twilight.  A  multicolor  photopola rime ter  scanned 
back  and  forth  in  azimuth  at  2.5  deg/sec,  incrementing  elevation  in  1* 
steps  between  5*  and  24*.  A  sample  of  these  data,  taken  in  March  1966 
at  5080A,  is  reduced  to  isolate  the  zodiacal  light.  Additional  data 
and  full  details  of  the  observations,  calibration,  and  data  reduction 
will  be  presented  elsewhere. 

Data  below  elevation  10*  are  omitted  here  due  to  difficulties  in 
the  atmospheric  corrections.  The  measured  brightnesses  were  converted 
to  absolute  units  (Sig(V)G„v)  by  reference  to  a  calibrated  (by  NBS  Fritz 
Peak  Observatory),  17.8-cm  diameter  1 ^-activated  phosphor  source.  The 
source  was  placed  over  the  objective  before  and  after  each  night's  ob¬ 
servations,  filling  both  the  aperture  and  the  3*  diameter  field  of  view 
(FOV).  Bright  stars  were  used  to  obtain  an  independent  absolute  cali¬ 
bration,  the  two  methods  agreeing  to  better  than  5  percent. 

Extinction  corrections  were  made  with  coefficients  derived  from  ob¬ 
servations  of  bright  stars  using  the  same  instrument.  Atmospheric  scat¬ 
tering  corrections  followed  the  method  outlined  by  Weinberg  (1964).  The 
brightness  contributed  by  "resolved"  stars  in  each  FOV  were  subtracted 
using  a  special  merged  star  catalog  developed  for  each  color  used  with 
this  instrument.  Background  starlight  was  subtracted  using  data  obtain¬ 
ed  from  Pioneer  10  observations  beyond  the  asteroid  belt,  whbre  the  godl- 
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ABSTRACT.  Background  starlight  observed  by  the  Pioneer  10  Imaging 
Photopolarimeter  from  beyond  the  asteroid  belt  la  used  to  Isolate 
zodiacal  light  in  Pioneer  observations  at  heliocentric  distances  B 
between  1  and  3  AU.  Near-ecliptic  zodiacal  light  brightness  data  in  the 
range  65°  to  l80°  elongation  e  are  used  to  depict  changes  In  the  shape 
of  the  zodiacal  light  with  e  and  R  and  are  compared  to  the  corresponding 
views  seen  from  the  Earth  and  from  the  Helios  1  and  2  spacecraft. 

The  heliocentric  dependence  of  zodiacal  light  can  only  be  observed 
directly  by  a  space  probe,  the  first  such  probe  measurement  being  made 
by  Pioneer  10.  During  the  cruise  phase  of  this  outer  solar  system  mis¬ 
sion,  sky  brightness  and  polarization  were  mapped  using  a  photopolari¬ 
meter  (Weinberg,  et  al.*,  1973,  197*0.  At  elongations  greater  than  90®, 
the  contribution  of  zodiacal  light  decreased  to  negligible  levels  beyond 
3  AU  (Hanner,  et  al.j  197k,  1976).  As  part  of  an  analysis  of  additional 
Pioneer  10  data,  Schuenuan,  et  al,  (1977)  pointed  out  that  there  is  no 
evidence ‘for  the  zodiacal  light  being  absent  beyond  3  AU  -  only  that  it 
becomes  vanishingly  small  compared  to  the  background  starlight  (integra¬ 
ted  starlight,  diffuse  galactic  light,  extragalactic  background  light). 
These  studies  are  extended  to  smaller  elongations  and,  for  the  first 
time,  zodiacal  light  brightnesses  are  tabulated  as  a  function  of  C.  8, 
and  R  (Table  l). 


Since  zodiacal  light  is  negligible  at  large  R,  observations  in 
these  regions  have  been  used  to  derive  maps  of  background  starlight  over 
the  sky  (Toller,  198li  Weinberg,  1981).  These  data  and  brightnesses  due 
to  discrete  (mv  <  6.9)  stars  are  subtracted  from  Pioneer  10  observations 
between  1  and  3  AU  to  derive  the  440oX  e.ft.R  zodiacal  light  topology 
given  in  Table  1.  The  three  entries  in  each  e,R  box  correspond  to  data 
at  ecliptic  latitudes  8  ■  +10°,  0°,  -10°  from  top  to  bottom,  respective¬ 
ly.  Empty  boxes  indicate  regions  not  yet  analyzed.  Blank  areas  within  a 
box  denots  regions  where  the  data  was  eitbsr  not  analyzed  or  is  unetr- 
taln  due  to  the  presence  of  a  bright  star  in  the  field  of  view. 
Additional  Pioneer  10  blue  data  and  the  corresponding  Pioneer  10  r«d 
(6400*)  data  will  be  similarly  analyzed  and  published  elsewhere. 
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ABSTRACT.  Recent  results  on  zodiacal  light  are  used  to  show  that  optics, 
dynamics,  and  infrared  must  be  considered  together  to  properly  and  fully 
characterize  the  interplanetary  dust  complex. 


Zodiacal  light  observations  have  been  widely  used  to  infer  the  large 
scale  properties  of  the  interplanetary  dust  -  with  mixed  results. 

Zodiacal  light  is  certainly  more  stable  (and  better  understood)  than  the 
published  literature  would  suggest  over  the  25  years  that  it  has  been 
studied  by  this  writer.  The  interested  reader  can  trace  this  history  in 
the  Proceedings  of  this  Colloquium's  predecessor  meetings  (Honolulu  1967, 
Heidelberg  1975,  Ottawa  1979),  in  triennial  Reports  on  Astronomy  of  IAU 
Commission  21,  and  in  various  reviews  (e.g.:  Leinert  1975 i  Weinberg  and 
Sparrow  1978;  Fechtig,  Leinert,  and  GrUn  1981).' 

There  have  been  intensive  ground  and  space  observations  over  the 
past  two  decades  together  with  laboratory  and  theoretical  studies,  and, 
as  we  entered  the  1980's,  there  was  general  agreement  that: 

-  zodiacal  light  has  solar  color  from  the  near  UV  to  the  near  IR,  except 
for  a  slight  reddening  in  shy  regions  observed  near  the  sun 

-  zodiacal  light  brightness  is  relatively  smooth  and  remarkably  stable 
over  times  as  long  as  a  solar  cycle  (Burnett  1976;  Dumont  and 
Levasseur-Regourd  1978;  Leinert,  et  al.  1982a) 

-  zodiacal  light  brightness  decreases  monotonically  with  heliocentric 
distance  R  and  is  negligible  beyond  the  asteroid  belt  (Hanner,  et  al.; 
1974,  1976) 

-  zodiacal  light  is  partially  plane  polarized  with  its  electric  vector 
perpendicular  to  the  scattering  plane,  except  for  regions  at  large 
elongations  where  there  is  polarization  reversal  (electric  vector 
parallel  to  scattering  plane). 

Recent  results  present  a  somewhat  different,  more  complex  picture 
of  the  interplanetary  dust.  Representative  of  this  changing  "view”  are 
the  recent,  exciting  observations  by  IRAS.  Figure  1  shows  IRAS  4-color 
data  on  total  thermal  emission  for  representative  ecliptic  pole-to-pole 
scans.  The  zodiacal  emission  peaks  near  the  ecliptic,  is  the  dominant 
source  at  12,  25,  and  60pm,  and  appears  to  be  the  brightest  diffuse 
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